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ABSTRACT 

O . 

We present a study of the spatial coherence of the intergalactic medium to- 
ward two pairs of high-redshift quasars with moderate angular separations ob- 
served with Keck/ESI, Q1422+2309A/Q1424+2255 (z em « 3.63, A9 = 39") and 
Q1439-0034A/B (z em « 4.25, AO = 33"). The crosscorrelation of transmitted 
flux in the Lya forest shows a 5 — la peak at zero velocity lag for both pairs. 
This strongly suggests that at least some of the absorbing structures span the 
230-300 /ifo 1 proper kpc transverse separation between sightlines. We also statis- 
tically examine the similarity between paired spectra as a function of transmitted 
flux, a measure which may be useful for comparison with numerical simulations. 

in 

teristics, we find that photon noise has little impact for S/N > 10 per resolution 
element. However, the agreement between the autocorrelation along the line 
sight and the crosscorrelation between sightlines, a potential test of cosmological 
geometry, depends significantly on instrumental resolution. Finally, we present 
an inventory of metal lines. These include a a pair of strong C IV systems at 
z ~ 3.4 appearing only toward Q1439B, and a Mg II + Fe II system present 

a' 



toward Q1439 A and B at z « 1.68. 

Subject headings: cosmology: observations — intergalactic medium — large scale 
structure of the universe - - quasars: absorption lines — quasars: individual 
(Q1422+2309, Q1424+2255, Q1439-0034A, Q1439-0034B) 
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1. 



Introduction 



Multiply-imaged lensed quasars and close quasar pairs provide valuable probes of struc- 
ture in the intergalactic medium. By comparing the absorption patterns in the spectra of 
adjacent quasar images one can gauge the similarity in the underlying matter distributions 
along the lines of sight and hence constrain the sizes of absorbing structures. Previous studies 
at small separations [AO < 10") have typically provided lower limits to the scale of Lya ab- 
sorbers over a wide range in redshift, with weak upper constraints of ~ 400 comoving kpc 
derived by assuming spherical clouds (Weymann & Foltz 1983; Foltz et al. 1984; McGill 1990; 
Smette et al. 1992; Dinshaw et al. 1994; Bechtold et al. 1994; Bechtold & Yee 1995; Smette et 
al. 1995; Fang et al. 1996). Observations of pairs at wider separations (A6 ps 0.5' — 3') have 
shown evidence that some Lyct absorbers span > 1 comoving Mpc (Petitjean et al. 1998; 
Crotts & Fang 1998; D'Odorico et al. 1998; D'Odorico, Petitjean, & Cristiani 2002; Young, 
Impey, & Foltz 2001; Aracil et al. 2002) and possibly up to 30 comoving Mpc (Williger et 
al. 2000). 

The majority of studies on lensed quasars and quasar pairs have relied on matching 
individual absorption lines between spectra. However, high- resolution spectra (FWHM < 20 
km s^ 1 ) are required to completely resolve these features. In addition, the crowding of 
lines in the high-redshift Lja forest often prevents the identification of single absorbers. An 
alternate approach reflecting the continuity of the underlying density field is to compute the 
correlation of transmitted flux along parallel lines of sight. This robust statistic is quickly 
computed and can be easily compared to numerical simulations of large-scale structure. 

As an extension to studying the matter distribution, correlations in the Lya forest 
have been proposed as a tool for constraining the cosmological constant through a variant 
of the Alcock-Paczyhski test (Alcock & Paczynski 1979). The Alcock-Paczyhski test takes 
advantage of the fact that, for a homogeneous sample of objects, the characteristic radial 
and transverse sizes should be equal. In the case of the Lya forest, one can compare the 
correlation length of absorbing structures along the line of sight to the correlation length 
in the transverse direction. The velocity separation, Av\\, between objects along the line of 
sight is simply given by their redshifts, 



In contrast, the transverse velocity separation, Av±, between objects at redshift z with angu- 
lar separation A9 depends on the cosmological parameters implicit in the Hubble constant, 
H(z), and angular diameter distance, D A (z), as 




Az 



(1) 



c. 



l + z 



Av ± = H(z)Al = H(z)D A (z)A6, 



(2) 
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where Al is the proper linear separation. One approach to exploiting the difference between 
Ai7|| and At>j_ is to directly compare the autocorrelation of transmitted flux along single lines 
of sight to the crosscorrelation between spectra of sources at a variety of angular separations 
(Hui, Stebbins, & Buries 1999; Lidz et al. 2003). An alternate method compares observed 
crosscorrelations to those determined from linear theory (McDonald & Miralda-Escude 1999) 
or from artificial spectra drawn from numerical simulations (Lin & Norman 2003). Recently 
Rollinde et al. (2003) found agreement between the crosscorrelations and autocorrelations 
among spectra of several z ~ 2 quasar pairs over a wide range in separation (although see 
the discussion on spectral resolution below). In the future, large surveys such as the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) should greatly increase the number of quasar 
pairs available for such studies. 

We present results for two pairs of quasars at a novel combination of high redshift 
and moderate separation. The closely separated A and C images (A9 = r','3) of the bright 
z = 3.63 lensed system QSO 1422+2309 (Patnaik et al. 1992) have been previously examined 
by Rauch, Sargent, & Barlow (1999, 2001) and Rauch et al. (2001). Adelberger et al. (2003) 
recently discovered an additional faint source, QSO 1424+2255 [z = 3.62) (therein referred 
to as Q1422b), at a separation A9 = 38'/5 from the lensed system. In this study, we compare 
the sightlines toward QSO 1422+2309A (herein referred to as Q1422) and QSO 1424+2255 
(herein referred to as Q1424 to avoid confusion with the B image of Q1422). We additionally 
investigate the A6 = 33'.'4 pair at z — 4.25 discovered in the Sloan Digital Sky Survey, SDSSp 
J143952.58- 003359.2 (herein Q1439A) and SDSSp J143951.60 - 003429.2 (herein Q1439B) 
(Schneider et al. 2000). For each pair, the similarity in quasar redshifts allows us to study the 
transverse properties of the Lya forest and intervening metal systems over a large pathlength. 
Our results should provide a valuable resource for comparison with numerical simulations of 
large-scale structure. 

The remainder of the paper is organized as follows: In §2 we present our observations 
together with a general overview of the data. We compute the flux correlation functions in 
the Lya forest in §3 and statistically examine the similarity between sightlines as a function 
of flux. In §4 we analyze the effects of photon noise and instrumental resolution on the 
correlation functions. Our results are summarized in §5, with an inventory of unpublished 
metal absorption systems presented in the appendix. 

Throughout this paper we adopt fl m = 0.3, = 0.7, and H = 70 km s _1 Mpc -1 . 
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2. The Data 

Our observations are summarized in Table 1. We observed all four quasars under good 
to excellent seeing conditions over the period 2000 March to 2002 June using the Keck 
Echellette Spectrograph and Imager (ESI) (Sheinis et al. 2002) in echellette mode. Additional 
observations of Q1424 were provided by C. Steidel while additional observations of Q1439A 
were provided by L. Hillenbrand. All exposures were taken at the parallactic angle except 
for one exposure of Q1439B at an airmass near 1.0, where chromatic atmospheric dispersion 
is only a minor concern. 

The raw CCD frames for Q1422, Q1439A, and Q1439B were processed and the 2-D 
echelle spectra extracted using the MAKEE software package. Reduction of Q1424 data was 
performed using a suite of IRAF scripts, as described in Adelberger et al. (2003). For each 
night we used the extracted orders from at least one standard star (Feige 34, BD+284211, 
and/or HZ 44) to determine an instrumental response function with which to derive relative 
flux calibrations. The calibrated orders of all exposures from all nights were then converted 
to vacuum heliocentric wavelengths and combined to produce a single continuous spectrum 
per object (Figures 1 and 2). We use a binned pixel size of 20 km s -1 . Normalized spectra 
were produced by fitting continua to the final reduced versions. 

The majority of observations were made using a 0V75 slit. However, inspection of the 
spectra extracted from exposures taken with different slits widths revealed very little dif- 
ference in resolution, likely due to favorable seeing. For the combined spectra we adopt 
a measured spectral resolution FWHM = 55 km s -1 . The final ESI spectra for Q1424, 
Q1439A, and Q1439B have typical S/N ^13 — 30 per resolution element in the Lya forest, 
while Q1422 has S/N w 150. 

We restrict our analysis of the Lya forest to the the wavelength region between each 
quasar's Lya and O VI emission lines. To avoid a proximity effect from the quasar, we include 
only those pixels at least 10,000 km s _1 blueward of the quasar's Lya emission. We further 
include only pixels at least 1,000 km s -1 redward of Ly/3 and O VI emission to avoid any 
confusion with the Ly/3 forest or intrinsic O VI absorption. For Q1422 and Q1424 this yields 
a redshift interval Az = 0.50 with mean redshift (z) = 3.22. For Q1439 A and B, Az = 0.58 
and (z) = 3.79. Due to the high redshift of Q1439 A and B, the Lya forest in the spectra 
of these sources extends over the strong night sky lines [O i] A5577 and Na I A5890,5896. In 
our analysis we exclude the narrow (2 — 3 A) regions around these lines. 

A visual comparison of the Lya forest in the paired spectra suggests that the most 
striking similarity occurs among the strongest absorption features (see Figures 3 and 4). 
These regions often appear to be similar along adjacent lines of sight, as do regions where 
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the absorption is nearly zero (possible "voids" relatively free of absorbing material). Matches 
among intermediate strength features are less obvious, which suggests that the structures 
giving rise to those lines may not be as coherent over the separation between sightlines. 
However, many strong lines and regions of nearly zero absorption that do not coincide 
between spectra can also be identified. The alignment of a subset of features may occur 
purely by chance. 

Our spectral coverage also allows us to investigate the extent to which metal systems 
span parallel lines of sight. Two strong C IV systems appear toward Q1439B at z ~ 3.4, 
separated by 1700 km s -1 . Each system has a large rest equivalent width, W rest (1548) ~ 1.5 
A. However, no C IV appears at this redshift toward Q1439A. In contrast, a low-ionization 
system containing Mg II and Fe II does appear along both sightlines at z ~ 1.68, separated by 
only ~ 400 km s -1 . A single, extended absorber responsible for the low-ionization features 
would have a linear size > 280 hj^ proper kpc. These lines may alternatively arise from the 
chance intersection of separate, clustered absorbers. Weaker C IV systems appear in both 
Q1422 and Q1424 separated by ~ 440 km s _1 at z ~ 3.08, where the separation between 
sightlines is 290 hj^ proper kpc. We expand further on the properties of metal systems in 
the appendix. 

In order to assess the impact of spectral characteristics on our Lya forest results we 
additionally employed a deep, high-resolution spectrum of Q1422+2309A. Observations were 
made with the Keck High Resolution Echelle Spectrometer (HIRES) (Vogt et al. 1994) using 
a / /574 slit. This yields a spectral resolution FWHM = 4.4 km s _1 . Reductions were 
performed as described Rauch, Sargent, & Barlow (2001). In addition to the data from that 
work, which used the red cross-disperser only, we include subsequent exposures taken with 
the UV-blazed cross-disperser installed to provide additional coverage of the Lya forest. The 
final spectrum was binned to give a constant velocity width for each pixel of 2.1 km s _1 , 
with a typical S/N per resolution element of 50 — 90. 

3. Comparison of Sightlines 

3.1. Correlation Functions 

The correlation of transmitted flux in the spectra of closely separated quasars provides a 
simple means of quantifying the degree of similarity in the matter distribution along adjacent 
sightlines. We define the unnormalized correlation, £, of the spectra of two sources separated 
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on the sky by an angle AO as 



£(A0, At;,,) = 1 J>(A0, W || + At;,,) - F(A0)][F(O, „„) - F(0)], 



(3) 



where F is the continuum-normalized flux, F(A9) and F(0) are the mean fluxes along the 
two sightlines, v\\ is the line-of-sight velocity, Av\\ is the longitudinal velocity lag, and N 
is the total number of pixels in each spectrum within the region of interest. The sum is 
performed over all available pixels at a given velocity lag, of which there will be 



for a pixel size Svp^. The normalized correlation value can be computed by dividing Eq. (3) 
by the standard deviation in each input spectrum. 

We have chosen a pixel size for the combined spectra to give roughly 3 pixels/resolution 
element. However, we note that the 1/N factor in Eq. (3) implies that the value of the 
correlation will be insensitive to pixel size so long as the spectrum is well sampled. The 
increase in the sum created by using a larger number of smaller pixels will be offset by the 
increase in N so long as there are > 2 pixels/resolution element. Using pixels larger than 
the spectral resolution will introduce smoothing effects (see §4). 

In order to assess the coherence of absorbing structures across adjacent sightlines, we 
compute the crosscorrelation of transmitted flux in the Lya forest in the spectra of our quasar 
pairs. As a reference, we also determine the autocorrelation, £(0, Av\\), along lines of sight 
toward individual objects. At these high redshifts, H I Lya absorption will strongly dominate 
over contaminating absorption from lower-redshift metals such as C IV and Mg II. The 
correlation functions should therefore accurately reflect the distribution of neutral hydrogen 
to within the present measurement errors. All correlations are computed in single-pixel steps, 
which are 20 km s" 1 for the ESI data. Figure 5 displays the autocorrelation function for 
each quasar along with the crosscorrelation functions between adjacent sightlines. Both in 
the case of Q1422/Q1424 and Q1439A/B, a clear peak in the crosscorrelation at zero lag 
indicates a genuine similarity between sightlines. 

Undulations in the correlation functions arising from the chance superposition of unre- 
lated lines constitute the dominant source of uncertainty in the peak values (see discussion 
on photon noise below). Pixel-to-pixel variations in the correlations are themselves clearly 
correlated. However, we find that the overall distributions of values away from the central 
peaks are very nearly Gaussian. We therefore take the standard deviation of pixels in the 
"noise" region, which we define to be where 2000 km s _1 < |Ai>||| < 18000 km s~ x , as the 
la error in a correlation peak value. In this region we expect no underlying signal, however 
the correlation is still computed from at least half of the available pixels. Since the sum 



n pix (Av\\) = N - Av\\/5v } 



pix 



(4) 
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in Eq. (3) is computed over fewer pixels as the velocity lag increases, yet the factor 1/N 
remains constant, the amplitude of the noise features will tend to diminish as ^n pix (Av\\) , 
where n p i x (Avn) is the number of pixels included in the sum in Eq. (3) at lag Av». Therefore, 
in order to match the amplitude of the noise at Av» = 0, we multiply the correlation at each 
lag by a scale factor s before computing the standard deviation, where 

/ rWO) _ / N 
\/n pix (At;||) y N-Av {l /Sv pix [) 

For the range in velocity lag defined above this is a modest correction. Including only lags 
where at least half of the ~ 1800 Lya pixels in each spectrum overlap limits s to at most 
\/2. For the range in velocity lag shown in Fig. 5, s < 1.2. 

Our results for the crosscorrelations are summarized in Table 2. Using the above es- 
timate for the error, the peak in the crosscorrelation for Q1422/Q1424 (Q1439A/B) is sig- 
nificant at the 7<j (5cr) level. This strongly suggests coherence in the absorbing structures 
on the scale of the 230 — 300 proper kpc transverse separation between sightlines. The 
marginal consistency of the peak values with one another likely reflects the similarity in 
sightline separation and redshift for the two quasar pairs. The zero-lag values of the au- 
tocorrelations, which give the variance in the flux for these sections of the Lya forest, are 
0.0950 for Q1422, 0.0796 for Q1424, 0.0862 for Q1439A, and 0.1000 for Q1439B. Thus, the 
normalized crosscorrelation peaks are 41.2 ± 5.7% for Q1422/Q1424 and 31.9 ± 6.2% for 
Q1439A/B. 



3.2. Flux Distributions 

The flux crosscorrelations demonstrate that at least some Lya absorbers span the sepa- 
ration between our paired lines of sight. However, it does not indicate whether the similarity 
between sightlines depends on the strength of the absorber. The crowded nature of the 
Lya forest at z ~ 4, together with the present spectral resolution, greatly inhibits a study of 
individual lines. However, we are able to look at the agreement between absorption features 
on a pixel-by-pixel basis. 

Our goal is to determine whether the similarity in flux between paired spectra depends 
on the amount of absorption for an individual pixel. First we first select those Lya forest 
pixels in one spectrum (Q1422 or Q1439A) whose flux falls within a specified range. We then 
identify the pixels in the companion spectrum with matching wavelengths and compute their 
flux distribution. Comparing the distribution in this subsample to that in all Lja forest pixels 
in the companion spectrum allows us to evaluate whether there exists an overabundance of 
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pixels in the specified flux range relative to that expected on random chance. 

The results for Q1422/Q1424 and Q1439A/B are shown in Figures 6 and 7, respectively. 
Each panel shows the normalized distribution of Q1424 or Q1439B pixels in the indicated 
subsample along with the distribution of all Lya forest pixels in that spectrum. The range 
of flux in Q1422 or Q 1439 A for each subsample is chosen to be significantly larger than 
the typical flux uncertainty. In each case, we compute the two-sided Kolmogorov-Smirnov 
statistic, which is the maximum difference between the cumulative fractions of pixels in the 
subsample and of all pixels in the forest, along with the associated likelihood of randomly 
obtaining a smaller statistic than the one observed. The number of pixels in each subsample 
is also shown. 

The clearest results occur at extreme levels of absorption. Pixels in Q1422 (Q1439A) 
either near saturation, F < 0.2, or near the continuum, F > 0.8, tend to strongly coincide 
with pixels of similar flux in Q1424 (Q1439B). Intermediate flux pixels appear to be less 
strictly matched, with the exception of pixels in Q1422 with 0.4 < F < 0.6. For certain 
subsamples, namely 0.6 < F < 0.8 in Q1422 and 0.4 < F < 0.6 in Q1439A, the distribution 
of fluxes in the companion spectrum is consistent with a random selection. This seems to 
confirm the visual appraisal that strong absorbers and regions relatively free of absorbing 
material tend to span adjacent lines of sight more readily than do absorbers of intermediate 
strength. 

While each sample includes enough pixels to produce a statistically significant result, a 
few caveats are worth considering. Systematic errors in continuum fitting might artificially 
create clusters of pixels with low absorption near the same wavelength in pairs of spectra. 
Similarly, errors in sky-subtraction that are repeated between spectra might create false co- 
incidences of pixels near zero flux. Multiple independent continuum fits resulted in typical 
differences of ~ 10% in flux on scales of ~ 100 A. We likewise find no evidence for large 
systematic errors in the sky subtraction. Given that the resulting uncertainties are signif- 
icantly smaller than the range in flux used to define a subsample in Figures 6 & 7, these 
effects should only be a minor concern. 

We stress that comparing transmitted fluxes does not strictly yield a clear physical in- 
terpretation. Pixels of intermediate flux commonly occur along the wings of strong features. 
They will therefore tend to cluster less readily than pixels near the continuum or near satu- 
ration, especially if the strong lines shift in velocity between spectra. Moreover, agreement 
in flux does not necessarily indicate agreement in optical depth. Since transmitted flux de- 
creases exponentially with optical depth, r, the difference in flux for a given fractional change 
in r will depend on the value of r itself. For a small characteristic change between sightlines, 
At/t <C 1, the greatest scatter in flux is expected for r ~ 1, with the scatter decreasing as 
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t — > or as t — > oo. Therefore, it is unclear whether the enhanced agreement in flux among 
pixels with flux near the continuum or near saturation indicates that the corresponding gas 
is more homogeneous on these scales than the gas giving rise to intermediate absorption 
features. More specific insights may be drawn by comparing our measured flux distributions 
to those derived from numerical simulations. 



4. Effects of Photon Noise and Instrumental Resolution 

Comparing the longitudinal and transverse flux correlation functions (or power spectra, 
equivalently) in the Lja forest has been been explored by several authors as a means of 
measuring the cosmo logical geometry via the Alcock-Paczyhski test. This application pri- 
marily requires a sample of sightlines large enough to overcome cosmic variance. However, 
some question remains regarding the dependence of the correlation functions on data char- 
acteristics such as signal-to-noise ratio and resolution. It is particularly important to know 
to what extent the autocorrelation measures the physical correlation length along the line 
of sight rather than the instrumental resolution. For additional discussion on the effects of 
spectral characteristics see, e.g., Lin & Norman (2003) and Croft et al. (1998). 

In the preceding analysis we assumed that the chance alignment of unrelated absorption 
features dominated over photon noise in producing uncertainty in the flux correlations. To 
justify this, we recomputed the autocorrelation for the ESI spectrum of Q1422 after adding 
increasing levels of Gaussian random noise. The results appear nearly identical for S/N 
> 10 per resolution element (Figure 8). A spike at Av\\ = appears in the autocorrelation 
function as the variance in the noise becomes comparable to the intrinsic variance in the 
absorption features. However, no such jump is expected to occur in the crosscorrelation 
since the photon noise in the two spectra should be uncorrelated. Large sets of moderate 
S/N spectra may therefore be more useful than smaller sets of high S/N for this type of 
work. 

To address the effects of spectral resolution on the autocorrelation we have synthesized 
moderate- and low-resolution spectra from a high-quality Keck HIRES spectrum of Q1422 
(resolution FWHM = 4.4 km s _1 ). In each test case we first tune the spectral resolution by 
convolving the HIRES data with a Gaussian kernel and then compute the resulting autocor- 
relation function. The results for spectra with resolution FWHM = 4.4 (unsmoothed), 15, 
55, and 200 km s _1 are plotted in Figure 9. Very little difference exists between the autocor- 
relation functions computed from the unsmoothed HIRES data and from the data smoothed 
to FWHM = 15 km s -1 . Thus, we may conclude that the correlation in the unsmoothed 
HIRES data is an accurate measure of the intrinsic correlation in the Lja forest (subject to 
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cosmic variance and redshift-space distortions). The Lya lines are already smoothed by their 
thermal width and easily resolved with HIRES. However, when the spectrum is degraded 
to FWHM = 55 km s -1 , similar to ESI data, the autocorrelation is clearly broadened and 
diminished in amplitude. (We note that the autocorrelation measured from this "synthetic" 
ESI spectrum is nearly identical to that computed from the real ESI data.) At even lower 
resolution, comparable to Keck/LRIS or VLT/FORS2, the spectral resolution dominates 
over the intrinsic correlation length. 

As the spectral resolution decreases, the peak of the autocorrelation also incorporates 
more of the outlying noise (random undulations in the correlation function at Av ^> 0). 
At large velocity lags (Avn ^> 200 km s -1 ), disagreement between the high- resolution and 
low-resolution cases may be due in part to these noise features. Better agreement may result 
when the correlation function is averaged over many sightlines. 

A more subtle issue is how spectral resolution affects the agreement between autocorre- 
lation and crosscorrelation functions. To address this, we degraded the resolution of our ESI 
data to mimic observations using a lower resolution instrument (FWHM = 200 km s _1 ) and 
then recomputed the correlations. Figure 10 contrasts the results from the original spectra 
with those from the smoothed versions. In each case we plot the autocorrelation functions 
computed from single sightlines and overplot the peak of the corresponding crosscorrela- 
tion at a velocity lag equal to the transverse velocity separation between the lines of sight 
(Av ± = 98 km s" 1 for Q1422/Q1422, Av ± = 96 km s' 1 for Q1439A/B for the case of Q m = 
0.3, Q A = 0.7, and H = 70 km s" 1 Mpc" 1 ). 

The differences shown in Figure 10 between the ESI and low-resolution cases suggest a 
significant dependence on spectral resolution. At ESI resolution the concordance between 
auto- and crosscorrelations appears to be very good for both pairs of sightlines. However, 
at low resolution the peaks of the crosscorrelation functions fall well below the values of the 
autocorrelation functions at the same total velocity separation. While part of this effect may 
be due to noise in the correlations, the general dependence on resolution can be understood in 
regions where the intrinsic correlation function is non-linear on scales smaller than the width 
of the smoothing kernel. (This includes the region around the peak at zero lag, since the 
correlation is expected to be symmetric about Avn = 0.) Convolving the input spectra with 
a smoothing kernel produces a correlation function that has been convolved twice with the 
same kernel, once for each spectrum. The resulting unnormalized autocorrelation function 
may therefore be higher or lower at a particular velocity lag, depending on the shape of the 
correlation function at that point. However, since the peak of the crosscorrelation function 
is already a maximum it can only decrease as a result of smoothing (apart from the effects of 
noise). Thus, if we adjust our cosmology such that, at high spectral resolution, the composite 
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crosscorrelation function built up from pairs at many different separations agrees with the 
mean autocorrelation function (ignoring redshift-space distortions), then this agreement may 
not hold when using low-resolution data. Spectral resolution must therefore be considered 
carefully when performing this type of comparison, for example by referring to artificial 
spectra drawn from numerical simulations of large-scale structure. 

5. Conclusions 

We have analyzed the transverse properties of the Lja forest at 3.0 < z < 4.1 on 
scales of ~ 1 h^Q comoving Mpc (230 - 300 /ifo 1 proper kpc) toward the z ~ 3.63 quasar 
pair Q1422+2309A/Q1424+2255 and the z « 4.25 pair Q 1439-0034 A/B. Strong peaks at 
zero velocity lag in the flux crosscorrelations between paired sightlines indicate coherence in 
the H I absorbers on these scales. The crowded nature of the Lya forest at these redshifts 
restricts a line-by-line comparison. However, a statistical approach suggests that the flux 
in paired spectra tends to be most similar in regions of either very strong or very weak 
absorption. The similarities at flux levels near the continuum are consistent with the large 
sizes of voids and low-density gaseous filaments and sheets expected to comprise much of 
the intergalactic medium. The more surprising similarities at low flux levels may simply 
reflect the limited dynamic range of transmitted flux in distinguishing among regions of high 
optical depth (r ^> 1). Larger differences in flux found at intermediate flux levels may relate 
to the increased scatter in flux naturally expected for optical depths near unity, and also to 
the fact that flux differences along the wings of strong lines will be sensitive to bulk velocity 
shifts between sightlines. Despite the difficulty in producing a clear physical interpretation 
from the data alone, both the correlation and flux statistics should provide a useful resource 
for comparison with numerical simulations. 

In order to assess the dependence of the correlation functions on photon noise and spec- 
tral resolution we have employed high-quality ESI and HIRES spectra of Q1422. Adding 
photon noise has very little effect for S/N > 10 per resolution element except to introduce a 
spike in the autocorrelation at zero lag. Increasing the resolution FWHM increases the veloc- 
ity width of the autocorrelation and decreases its amplitude, as expected. These effects may 
occur unevenly, however, due to the incorporation of noise features as the peak broadens. 
We further tested the effects of resolution on the relative amplitudes of the autocorrela- 
tion and crosscorrelation by degrading the resolution of our ESI spectra to match that of a 
low-resolution spectrograph such as Keck/LRIS or VLT/FORS2. For both pairs, the cross- 
correlation and autocorrelations agree well for our adopted cosmology when the correlations 
are computed from the original ESI data. However, smoothing the spectra substantially 
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reduces the crosscorrelation with respect to the autocorrelations. This will typically occur 
at any velocity lag where the intrinsic correlation function is non-linear on scales smaller 
than the width of the smoothing kernel. Cosmological tests making use of the Lya forest 
correlation length should take this effect into account. 

Finally, we present an inventory of metal lines in Q1424, Q1439A, and Q1439B in the 
appendix. The most noteworthy features are a pair of strong C IV complexes at z ~ 3.4 
appearing only toward Q1439B. Given their large equivalent widths, the relatively small sep- 
aration between these systems suggests that they are related, possibly arising from separate 
galaxy halos or two sides of a galactic outflow. Low-ionization systems containing Mg II and 
Fe II appear at z rs 1.68 in Q1439 A and B separated by 410 km s _1 . These may either 
represent a single, extended structure (A/ > 280 proper kpc) or separate, clustered ab- 
sorbers. A C IV system appears toward Q1424 at z — 3.084, within 440 km s" 1 of a known 
C IV system toward Q1422. A single absorber responsible for both systems would have a 
transverse size > 290 hjQ proper kpc. However, the number density of weak C IV systems 
at this redshift suggests that these lines may instead result from chance superposition. 
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Lynne Hillenbrand for kindly sharing her observations of Q 1439 A, and Rob Simcoe for his 
HIRES reduction of Q 1422 A. We would also like to thank Tom Barlow for making public 
the MAKEE software package, Bob Carswell for VPFIT, and the anonymous referee for 
helpful comments. Finally, we are indebted to the indigenous Hawaiian community for the 
opportunity to conduct observations from their sacred mountain. Without their hospitality 
this work would not have been possible. W.L.W.S gratefully acknowledges support from the 
NSF through grants AST-9900733 and AST-0206067. M.R. has been supported by the NSF 
through grant AST-0098492 and by NASA through grant AR-90213.01A. 

A. Metal Systems 

Previous studies using quasar pairs and multiply-imaged lensed quasars have demon- 
strated spatial coherence in metal systems over a variety of length scales. Differences between 
metal systems in narrowly-separated lines of sight suggest that the absorbers giving rise to 
individual lines seen in high-resolution spectra span at most a few kiloparsecs, both for high- 
ionization systems seen in C IV and for low-ionization systems seen in Mg II (Lopez et al. 
1999; Petitjean et al. 2000; Rauch, Sargent, & Barlow 2001; Rauch et al. 2002; Churchill et 
al. 2003; Tzanavaris & Carswell 2003; Ellison et al. 2004). However, these absorbers may 
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be part of larger structures extending > 20 1i 7 q kpc (Smette et al. 1995) and even > 100 
h^Q kpc for highly ionized material (Petitjean et al. 1998; Lopez et al. 1999; Lopez, Hagen, 
& Reimers 2000). 

The present data afford us a unique opportunity to probe the coherence of metal systems 
on scales of ~ 300 proper kpc. Unfortunately, modest signal-to-noise and significant 
contamination from skylines in the red part of our spectra greatly limit our sensitivity and 
hinder completeness estimates. Our detections are limited to relatively strong lines and lines 
that fall in regions of unusually high S/N. 

Line lists for Q1424, Q1439A, and Q1439B are presented in Tables 3, 4, and 5, respec- 
tively. In the following sections we briefly comment on some of the more interesting systems. 
Results for Q1422 are presented in detail elsewhere (Rauch, Sargent, & Barlow 2001; Bech- 
told & Yee 1995; Petry, Impey, & Foltz 1998). For each line we measure the wavelength 
centroid and equivalent width. Weighted mean redshifts are given where multiple lines are 
measured for a single ion. In the case of blended lines we attempt to alleviate the overlap by 
fitting Voigt profiles to any unblended transitions of the same ions using VPFIT and then 
dividing by the inferred model profiles. This typically allows us to obtain values for at least 
the strongest blended components. 



A.l. Q1424+2255 

Strong associated broad absorption lines (BALs) at z — 3.62 are seen in C IV and 
N v along with weaker associated absorption in Si IV (although the Si IV appears affected by 
skyline contamination). Isolated components of C IV and N V appear up to 1000 km s" 1 blue- 
ward of the BAL. These suggest small pockets of outflowing material moving at a higher 
velocity than the material responsible for the bulk of the broad absorption. C IV is also 
seen at z — 3.62 toward Q1422A. (Rauch, Sargent, & Barlow 2001). However, since the 
absorption is associated with the quasar in both cases, these features are unlikely to be 
related. 

An intervening C IV system appears at z — 3.084. We additionally identify the line at 
5692.3 A to be Si iv A1394 at the same redshift, where Si iv A1403 is blended with N v A1239 
at z — 3.621. A C IV complex with possible Si IV appears at z — 3.090 toward both the 
A and C images of Q1422 (Rauch, Sargent, & Barlow 2001; Boksenberg, Sargent, & Rauch 
2003), implying a transverse size > 290 hjQ proper kpc for a single structure spanning all 
three lines of sight. Figure 11 shows the C IV absorption at this redshift in Q1424 and 
Q1422. The systems are narrowly separated along the line of sight, with Av rs 440 km s -1 . 
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However, Boksenberg, Sargent, & Rauch (2003) find ~ 25 C IV systems per unit redshift at 
z = 3 for systems with total column density N syst (C iv) > 10 12 cm -2 . This translates to a 
typical spacing between systems of only ~ 3000 km s -1 . It is therefore reasonable to suspect 
that this superposition of relatively modest C IV absorbers toward Q1422 and Q1424 occurs 
by chance. Significant H I absorption occurs along both lines of sight in the velocity interval 
spanned by these systems, though this may be due to the crowded state of the Lja forest 
at z ~ 3. 



A. 2. Q1439-0034A & B 

Two broad C IV complexes at z ~ 3.4 form the most conspicuous metal absorption 
features in the spectrum of Q1439B (Figure 12). The systems are separated by 1700 km s _1 , 
and each has a velocity width of ~ 400 km s -1 . Neither complex is observed toward Q1439A. 
The individual components are unlikely to be resolved with ESI. However, we can place a 
lower limit on the column density of these systems by assuming that the weaker transition 
(C IV A1551) falls on the linear part of the curve-of-growth. In that case, the column density, 
N syst , is related to rest equivalent with, W TCSt , as 

7T6 2 

West = ^/A 2 iV syst , (Al) 

m e c 2 

for a transition with rest wavelength Ao and oscillator strength /. For the lower-equivalent 
width system (z abs = 3.402, W res t(1551) « 1.1 A), this gives A" syst (C iv) > 5 x 10 14 cm" 2 . 
If we adopt a power-law column density distribution d 2 J\f / dN syst dz oc N~^ t , where M is the 
number of systems, then we would expect < 0.6 systems per unit redshift in this column 
density range based on the above Boksenberg et al. (2003) results and adopting their value of 
f3 — 1.6. The narrow redshift spacing of these two systems (Az = 0.025) therefore strongly 
suggests that they are related. One explanation is that the absorption complexes arise from 
two sides of an intervening galactic outflow. Alternatively, the high levels of H I absorption 
seen at both C IV redshifts may indicate a pair of galaxy halos. 

Mg II + Fe II systems appear toward both Q1439 A and B at z ~ 1.68, separated between 
spectra by 410 km s" 1 along the line of sight (Figure 13). From Churchill et al. (1999), the 
number of Mg II systems per unit redshift with rest equivalent width W / res t(2796) > 0.02 
A is expected to be ~ 2.7. The probability of randomly finding a system within this velocity 
separation is therefore ~ 0.02, or a factor of three lower if we restrict ourselves to strong 
(W rest (2796) > 0.3) systems (Steidel & Sargent 1992). Thus, there is a high likelihood that 
these systems are related. Both Mg II and Fe II appear stronger toward B than toward A, 
with a factor of seven difference between sightlines in the equivalent widths of Mg II A2804 
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and Fe II A2384. These features may arise either from a single absorber spanning the 280 
hjQ proper kpc between sightlines or from a pair of clustered objects. Two additional strong 
Mg II systems appear along only a single sightline, one toward A at z = 2.266 and another 
toward B at z = 1.399. 

Q1439B is a BAL quasar showing high levels of associated absorption in C IV and N v. 
We identify the O I absorption reported by Fukugita et al. (2004) as the strong C IV systems 
at z ~ 3.4 mentioned above. No strong associated absorption is seen for Q1439A. As 
noted by Fukugita et al. (2004), the difference in associated absorption features supports the 
hypothesis that Q1439 A and B are a true binary and not a lensed pair. 
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Table 1. Summary of Observations 



Object Date a Exp. time FWHM b 

(s) (km s- 1 ) 



Keck/ESI 



Q1422+2309A 


2000 Mar 03-04 


1200 c 


~ 55 


Q1424+2255 


2000 Mar 03 - 2001 Apr 21 


44900 d 


~ 55 


Q1439-0034A 


2001 Jan 26 - 2002 Jun 10 


8000° 


~ 55 


Q1439-0034B 


2001 Apr 19 - 2002 Jun 10 


42400 f 


~ 55 


Keck/HIRES 


Q1422+2309A 


1998 Jan 30 - 1998 Apr 15 


31600 s 


4.4 



a Interval over which observations were made. 
b Spectral resolution. 

c 600 s with ff.'75 slit, 600 s with l'.'O slit. 

d 43100 s with 0C75 slit (including 20900 s by C. Steidel), 1800 s with 
l'.'O slit. 

e 5000 s with 0'.'75 slit, 3000 s with (X'5 slit (L. Hillenbrand). 
f 0"75 slit. 

g 9000 s with UV cross-dispcrser, 22600 s with red cross-disperser 



Table 2. Lya Forest Crosscorrelation Values 



QSO Pair 


A6» a 


z Lya h A/ c 


Av ± d 


«A0,O) 






h^Q proper kpc 


km s _1 


Unnormalizcd 


Normalized* 


Q1422/Q1424 
Q1439A/B 


38'.' 5 
33'.' 4 


2.965 - 3.466 283 - 298 
3.496 - 4.075 230 - 245 


97.8 
96.3 


0.0358 ± 0.0049 
0.0296 ± 0.0057 


0.412 ±0.057 
0.319 ±0.062 



a Angular separation of the QSO pair. 

b Rcdshift interval used in computing the correlation values for the Lya forest. 

c inear separation between the lines of sight for the given redshift interval, computed for Q m = 0.3, 
Q A = 0.7, and H = 70 km s" 1 Mpc" 1 . 

d Transverse velocity separation between the lines of sight for our adopted cosmology at the mean redshift 
in the given interval. 

e Crosscorrelation at zero lag computed from Eq. (3). 

'Crosscorrelation at zero lag computed from Eq. (3) divided by the standard deviation of the Lya forest 
flux in each input spectrum. 
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Table 3. Q 1424+2255 Metal Systems 



A obs Line ID z obs a W rest 

(A) (A) 



Intervening 



5692.34 ±0.12 


Si iv A1394 


3.08417 ±0.00009 


0.29 ± 0.02 


6323.34 ±0.14 


C iv A1548 


3.08433 ± 0.00007 


0.50 ± 0.03 


6333.93 ±0.17 


C iv A1550 


3.08433 ± 0.00007 


0. 


37 ±0.03 


Associated 


5704.76 ± 0.08 


N v A1239 


3.60499 ± 0.00007 


0. 


34 ±0.01 


5709.70 ±0.11 


N v A1239 


3.60898 ± 0.00009 


0. 


11 ± 0.01 


5724.81 ± 0.08 


N v A1239 


3.62118 ±0.00005 


0. 


96 ± 0.02 


5743.22 ±0.09 


N v A1243 


3.62118 ±0.00005 





75 ± 0.02 


6437.34 ±0.17 


Si iv A1394 


3.61865 ± 0.00009 





24 ± 0.03 


6478.79 ±0.21 


Si iv A1403 


3.61865 ± 0.00009 


0. 


16 ± 0.02 


7129.64 ±0.08 


C iv A1548 


3.60511 ±0.00005 b 





36 ± 0.01 


7135.69 ±0.16 


C iv A1548 


3.60902 ± 0.00010 


0. 


10 ±0.01 


7141.28 ±0.08 


C iv A1550 


3.60511 ±0.00005 b 




< 0.37 c 


7144.70 ±0.13 


C iv A1548 


3.61484 ± 0.00008 





08 ± 0.01 


7153.99 ±0.07 


C iv A1548 


3.62085 ± 0.00004 


1 


82 ± 0.03 


7165.96 ±0.11 


C iv A1550 


3.62085 ± 0.00004 


1. 


44 ± 0.06 



a Weighted mean redshift from all available transitions. 
b C IV doublet redshift computed from the A1548 line only. 
c 3cr upper limit (possible blend). 



Table 4. Q1439-0034A Metal Systems 



A bs Line ID z obs a W res t 

(A) (A) 



Intervening 



6393.06 ±0.10 Fc II A2383 1.68304 ± 0.00004 0.07 ± 0.01 

7502.82 ±0.12 Mg n A2796 1.68305 ± 0.00004 0.18 ±0.02 

7521.87 ±0.19 Mg II A2804 1.68305 ± 0.00004 0.10 ±0.02 

9134.29 ±0.14 Mg II A2796 2.26645 ± 0.00004 0.30 ± 0.03 

9157.37 ±0.18 Mg II A2804 2.26645 ± 0.00004 0.32 ± 0.04 



a Weighted mean redshift from all available transitions. 
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Table 5. Q1439-0034B Metal Systems 



A bs Line ID z obs a W res t 

(A) (A) 



Intervening 



6402 


.03 


± 


0. 


06 


Fe II A2382 


1 


.68680 


± 





.00002 





.51 


± 





.01 


6602. 


.02 


± 


0. 


13 


C iv A1548 


3 


.26419 


± 





.00006 





15 


± 





.02 


6612 


.57 


± 


0. 


14 


C iv A1551 


3 


.26419 


± 





.00006 





.13 


± 





.02 


6704. 


.31 


± 


0. 


11 


C iv A1548 


3 


.33041 


± 





.00006 





.14 


± 





02 


6707. 


.40 


± 


0. 


06 


Mg II A2796 


1 


.39863 


± 





.00002 


1 


.18 


± 


0. 


.03 


6709. 


.91 


± 





11 


Si II A1527 


3 


.39502 


± 





.00007 





.17 


± 





.02 


6715 


.59 


± 


0. 


19 


C iv A1551 


3 


.33041 


± 





.00006 





.09 


± 





.02 


6724. 


.68 


± 


0. 


07 


Mg II A2804 


1 


.39863 


± 





.00002 





.91 


± 


0. 


.03 


6804. 


02 


± 


0. 


07 


C iv A1548 


3 


.39476 


± 





.00003 


1 


.48 


± 





02 


6815. 


.18 


± 


0. 


08 


C iv A1551 


3 


.39476 


± 





.00003 


1 


.09 


± 


0. 


.02 


6843 


.29 


± 


0. 


06 


C iv A1548 


3 


.42022 


± 





.00003 


1 


.58 


± 





.04 


6854 


.89 


± 


0. 


06 


C iv A1551 


3 


.42022 


± 





.00003 


1 


.25 


± 





.02 


6949. 


.75 


± 


0. 


20 


Fe ii A2587 


1 


.68680 


± 





.00002 





.26 


± 





.05 


6986. 


.15 


± 





08 


Fe ii A2600 


1 


.68680 


± 





.00002 





.50 


± 





02 


7387 


.85 


± 





14 


C iv A1548 


3 


.77191 


± 





.00007 





18 


± 





02 


7400. 


22 


± 





19 


C iv A1551 


3 


.77191 


± 





.00007 





.09 


± 





02 


7513 


.04 


± 


0. 


08 


Mg ii A2796 


1 


.68673 


± 





,00003 b 





.70 


± 





.03 


7532. 


.51 


± 





15 


Mg ii A2804 


1 


.68673 


± 





.00003 b 





.68 


± 





.13 



Associated 



6502.00 ± 0.09 


N v A1239 


4.24877 ± 0.00005 


1.62 ±0.05 


6523.45 ± 0.08 


N v A1243 


4.24877 ± 0.00005 


1.31 ±0.02 


8125.49 ± 0.04 


C iv A1548 


4.24847 ± 0.00002 


1.94 ±0.01 


8139.44 ± 0.04 


C iv A1551 


4.24847 ± 0.00002 


1.68 ± 0.01 



a Wcighted mean rcdshift from all available transitions. 
b Mg II doublet redshift computed from the A2796 line only. 
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5000 5500 6000 6500 7000 



A (A) 

Fig. 1. — Keck/ESI spectra covering the rest wavelength region from Ly(3 to C IV in Q1422 
(top) and Q1424 (bottom). Spectra have been binned in wavelength using 40 km s _1 pixels 
for display. Contamination from the [O i] A5577 skyline can be seen in the spectrum of 
Q1424. 
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A (A) 

Fig. 2.— Same as Figure 2, here for Q1439A (top) and Q1439B (bottom). 
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Fig. 3. — A representative section of the Lya forest in the normalized spectra of Q1422 (grey- 
line) and Q1424 (black line). 
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Fig. 5. — Flux correlation functions computed from the Lyct forest plotted vs. longitudinal 
velocity lag. Top four panels: Autocorrelation functions for Q1422 (top left), Q1424 (middle 
left), Q1439A (top right), and Q1439B (middle right). Negative velocity lags have been in- 
cluded for consistency. Bottom panels: Crosscorrelation functions for Q1422/Q1422 (bottom 
left) and Q1439A/B (bottom right). 
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Q1424 



Fig. 6. — Distributions of continuum-normalized pixel fluxes in the Lya region of Q1424 as 
a function of the flux in Q1422. Unshaded histograms show the flux probability distribution 
for the subsample of pixels in Q1424 corresponding in wavelength to those pixels in Q1422 
that have fluxes in the indicated range. Shaded histograms show the probability distribution 
for all Q1424 pixels in the Lya region. For each range of flux in Q1422, the two-sided 
Kolmogorov-Smirnov statistic for the subsample and full sample of Q1424 flux values is 
presented along with the associated likelihood, expressed as a percent, of obtaining a higher 
value if the the two samples were drawn from the same distribution. The number of pixels 
in each subsample is given in parentheses. The full sample contains 1783 pixels. 
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Fig. 7. — Same as Figure 6, here for pixel fluxes in the Lya region of Q1439B as a function 
of the flux in Q1439A. The full sample contains 1793 pixels. 
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200 400 600 

AVm (km/s) 

Fig. 8. — The effects of photon noise on the shape of the autocorrelation function. The solid 
line shows the autocorrelation computed from our original ESI spectrum of Q1422 (FWHM 
= 55 km s _1 , S/N = 150 per resolution element). Additional lines show the autocorrela- 
tion recomputed after adding random Gaussian noise to the original spectrum such that 
the resulting S/N per resolution element is 20 (dash-dotted line), 10 (dashed line), and 5 
(dotted line). Autocorrelation functions were computed in longitudinal velocity lag steps of 
20 km s _1 . The spike at Av» = reflects only a single point for each curve. 
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Fig. 9. — The effects of spectral resolution on the shape of the autocorrelation function. 
The solid line shows the autocorrelation computed from our original HIRES spectrum of 
Q1422 (FWHM = 4.4 km s^ 1 ). Additional lines show the autocorrelation recomputed after 
smoothing the HIRES data to a spectral resolution FWHM of 15 km s _1 (dash-dotted line), 
55 km s _1 (dashed line), and 200 km s _1 (dotted line). 
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Fig. 10. — The relative effects of spectral resolution on auto- and crosscorrelations. Top pan- 
els: Autocorrelation functions computed from the Lya regions in the indicated unsmoothed 
ESI spectra (FWHM = 55 km s -1 ) (solid and dotted lines) overlaid with the peak value of 
the crosscorrelation computed between those spectra (filled circles). The peak of the cross- 
correlation is plotted at a velocity lag equal to the transverse separation between the two 
lines of sight for the case of fl m = 0.3, ^a= 0.7, and H = 70 km s" 1 Mpc -1 . Bottom panels: 
Same as the top panels, here after smoothing the ESI spectra to a resolution of 200 km s _1 . 
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Fig. 11. — Top three panels: Continuum-normalized ESI spectrum of Q1422 (histogram) 
showing the C IV and H I absorption for the system at z = 3.0902 reported by Rauch, Sargent, 
& Barlow (2001). The continuous lines show the Voigt profile fits computed from HIRES data 
where the profiles have been smoothed to ESI resolution. Bottom three panels: Continuum- 
normalized ESI spectrum of Q1424 (histogram) showing the C IV and H I absorption for the 
system measured herein at z = 3.0843. The continuous lines show Voigt profile fits computed 
from the ESI data. All components in both spectra are placed on a common velocity scale 
in order to demonstrate the longitudinal velocity separation between these systems. 
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Fig. 12. — Top three panels: Continuum-normalized ESI spectrum of Q1439A (histogram) 
covering the C IV and H I regions near z = 3.4. The continuous lines in the top two 
panels demonstrate the absence of C IV over this interval. Bottom three panels: Continuum- 
normalized spectrum of Q1439B (histogram) showing the C IV and H I absorption for the 
systems at z = 3.3948 and z = 3.4202. The continuous lines show the Voigt profile fits to the 
C IV absorption. Several of the C IV A1548 pixels near 850 km s _1 in this plot are heavily 
affected by skylines and have been excluded from the fit. All components in both spectra 
are placed on a common velocity scale in order to demonstrate the longitudinal velocity 
separation between these systems. 
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Fig. 13. — Top three panels: Continuum-normalized spectrum of Q1439A (histogram) show- 
ing the Mg II and Fe II absorption for the system at z — 1.6831. The continuous lines 
show Voigt profile fits. Bottom three panels: Same as the top panels, here for the system at 
z = 1.6867 in Q1439B. Several of the Mg II A2804 pixels in Q1439B near 200 km s" 1 in this 
plot are heavily affected by skylines and have been excluded from the fit. All components in 
both spectra are placed on a common velocity scale in order to demonstrate the longitudinal 
velocity separation between these systems. 



